function F = mply_prft(prc,nu,mc,a,b,par_beta)

%
% This program computes the monopoly profit given prices
% Code relies on the fact that outside good is the first element
% of the price, nu and mc vectors.  Price is set of zero, nu to 
% the choosen parameter and mc =0;

prc_guess = [0;prc];

% demand
     dmd = demand(prc_guess,nu,a,b,par_beta);
     
% profit
     pi_vec = (prc_guess - mc).*dmd;
     
%     [prc_guess dmd pi_vec], pause
     
     F = -sum(pi_vec);

